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Differential elliptic flow cjHfrito pj_ ~ 5 GeV/c is calculated using the MPC elas- 
tic parton cascade modcoEl for Au+Au at E cm ry 130 A GeV. The results are 
compared to recent STAR/RHIC elliptic flow datatl. An elliptic flow pattern com- 
parable to the data seems to require initial parton densities at least twice higher 
than the predicted by the HIJING modelta. Elliptic flow is also shown to be sensi- 
tive to the hadronization procedure. 



1 Introduction 

Elliptic flow, V2(p±) = (cos(2(j)))p ± , the differential second moment of the 
hal momentum distribution, has been the subject of increasing interest 
since the discovery!] at RH1C that V2(p±_ > 1 GeV) — > 0.2. This sizable 
high p± collective effect depends strongly on the dynamics in a heavy ion 
collision and provides important information about the density and effective 
energy loss of partons. 

The simplest theoretical framework to study elliptic flow is ideal 
hydrodynamicsa. For RHIC energies, ideal hydrodynamics agrees remark- 
ably well with the measured elliptic flow dataB up to transverse momenta 
~ 1.5 GeV/c. However, it fails to saturate at high p± > 2 GeV as does the 
data reported by STAR at Quark Matter 2001. 

A theoretical problem with ideal hydrodynamics is that it assumes local 
equilibrium throughout the whole evolution. This idealization is marginal 
for conditions encountered in heavy ion collisionsB. A theoretical framework 
is required that allows for noncquilibrium dynamics. Covariant Boltzmann 
transport theory provides a convenient framework that depends on the local 
mean free path X(x) = l/an(x). It interpolates between free streaming (A = 
oo) and ideal hydrodynamics (A.gd}). 

Parton cascade simulations^' 3tL3 show on the other hand, that the initial 
parton density based on HIJINGa is too low to produce the observed elliptic 
flow unless the pQCD cross sections are artificially enhanced by a factor ~ 

2 — 3. However, gluon saturation modelstll predict up to five times higher 
initial densities, and these may be dense enough to generate the observed 
collective flow even with pQCD elastic cross sections. In this study, we explore 
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the dependence of elliptic flow on the initial density, or equivalently^j, on the 
partonic cross section. __ 

Calculations based on inelastic parton energy lossacl also predict satura- 
tion or decreasing V2 at high p± . These calculations are only valid for high p± , 
where collective transverse flow from lower p± partons can be neglected. The 
hydrodynamic component from low p± is, on the other hand, automatically 
incorporated in parton cascades. Though parton cascades lack at present co- 
variant inelastic energy loss, we find in Fig. [I] that elastic energy loss alone 
may account for the observed high p± azimuthal flow pattern as long as the 
elastic opacity is large enough. _ 

We compute the partonic evolution with MPCo, a newly formulated, co- 
variant, parton kinetic theory technique. MPC is an extension of the covariant 
parton-cascade algorithm, ZPGL3. Both MPC and ZPC have been extensively 
testect3o and compared to analytic transport solutions and covariant Euler 
and Navier-Stokes dynamics in 1+1D geometry. A critical new element ofjjipjjh 
these algorithms is the parton subdivision technique proposed by PangE 

2 Covariant Parton Transport Theory 



We consider here, as in Refs.ul3E3lIj, the simplest form of Lorentz-covariant 
Boltzmann transport theory in which the on-shell phase space density f(x,p), 
evolves with an elastic 2 — > 2 rate as 

P^h = jjjihh ~ hh)W 12 ^ A 8\p 1 + P2 - P3 -p 4 ) + S(x, Pl ). (1) 

234 

Here W is the square of the scattering matrix element, the integrals are short- 
hands for J = J MnWlt- ' wnere 9 is the number of internal degrees of freedom, 

i 

while fj = f(x,pj), and we interpret /(ir,p) as describing an ultrarelativistic 
massless gluon gas with g = 16 (8 colors, 2 helicities). The initial conditions 
are specified by the source function S(x, p), which we discuss in Section ^. 

The elastic gluon scattering matrix elements in dense parton systems are 
modeled by a Debye-screened form 

dt V s J (t - m 2 ) 

where [i is the screening mass, Co = 97ra 2 /2/i 2 is the total cross section, which 
we chose to be independent of energy. For a fixed <7o, the relevant transport 
cross section cr t — J dersin 2 9 cm is maximal in the isotropic (/x — > oo) case. 

a The equivalence is due to the scaling property explained in Section [j| 
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It is important to emphasize that while the cross section suggests a geo- 
metrical picture of action over finite distances, Eq. (|^) is only a convenient 
parametrization to describe the effective local transition probability, W , as 
dW/dt — s dajdt. The particle subdivision technique (see next Section) re- 
moves all notion of nonlocality in this approach, as in hydrodynamics. 

3 Parton Subdivision and Scaling of Solutions 

Cascade algorithms inevitably violate Lorentz covariance because particle in- 
teractions are assumed to occur whenever the distance of closest approach (in 
the relative cm.) is d < \J o-q/h, which corresponds to action at a distance. 
To recover the local character of equation. Inf) and hence Lorentz covariance, 
we use the parton subdivision techniquecj'Ej. This technique is based on the 
covariance of Eq. ([l]) under the transformation 

/' = //, W^W' = W/l (a - a' = tr/Z), (3) 

which effectively rescales the cascade interaction range by Lorentz 
violation formally vanishestil in the_J — * oo limit. In practice, very high 
~ 100 — 1000 subdivisions are neededo to obtain accurate numerical solutions 
from initial conditions expected at RHIC. 

Subdivision covariance (^|) implies that the transport equation has a broad 
dynamical range, and the solution for any given initial condition and trans- 
port property immediately provides the solution to a broad band of suitably 
scaled initial conditions and transport properties. This is because solutions 
for problems with / times the initial density dN/drj, but with one l-th the 
reaction rate can be mapped to the original (I — 1) case. We must use sub- 
division to eliminate numerical artifacts. However, once that is achieved, we 
have actually found the solution to a whole class of suitably rescaled problems. 

The dynamical range of Eq. (|l|) is further increased by its covariance 
under coordinate and momentum rescalingcl, leading to covariance under 

^(fe})^ W '({M) = ^^({^}) 5 m^ m i = m/l p , (4) 

where l x and l p arc the coordinate and momentum scaling parameters, respec- 
tively. This meansB that we can scale one solution to others provided that 
/i/Xb, aodN/drj and Rq/tq remain the same. 
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4 Results and Discussion 
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Figure 1. The proper time evolution of elliptic flow as a function of transverse momentum, 
for an Au+Au collision at t/s = 130A GeV with b = 8 fm, in case of isotropic (top row) 
and gluonic (bottom row) cross section with no = 3 mb (left) and 10 mb (right). 

Following Ref.i, we modeled the colliding nuclei as longitudinally boost 
invariant Bjorken cylinders with radii Rq — 6 fm in local thermal equilib- 
rium at temperature To at proper time To = 0.1 fm/c. The pseudo-rapidity 
i] = 1/2 log((i + z)/(t— z)) distribution was taken as uniform between \r/\ < 5. 
For central collisions, T = 500 MeV and dN/drj = 400 as by fitting the 
gluon mini-jet transverse momentum spectrum predicted by HIJINGu (in- 
cluding shadowing and jet quenching effects). For non-central collisions, as a 
crude estimate, we scaled dN/drj by the overlap area of the two cylinders. To 
maximize elliptic flow, the impact parameter was chosen to be b = 8 fm. 

Evolutions from different initial densities can be obtained by varying the 
cross section only and using the scaling property explained in Section |^. Our 
simulations were performed both with isotropic and gluonic cross sections 
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(/j = 0.5 GeV) with <tq = 3 and 10 mb. The particle subdivision was I = 100. 

Fig. |l| shows the proper time evolution of the gluon elliptic flow as a 
function of transverse momentum for the different cross sections. For all cross 
sections studied, elliptic flow reaches its final asymptotic value early, by r = 4 
fm/c, reinforcing results in Ref. □. Also, elliptic flow increases with increasing 
transport cross section, as expected. 

Elliptic flow monotonically increases until p± ~ 2 GeV, where it satu- 
rates. According to a more detailed studylilJ, this saturation occurs for all 
impact parameters, leading to a saturation in the impact parameter averaged 
(minimum bias) elliptic flow as well. For the 3 mb gluonic cross section, t>2 
decreases after 2 GeV, due to elastic energy loss. 
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Figure 2. Partonic and the final hadronic elliptic flow after hadronization via independent 
or Lund fragmentation, for the 10 mb gluonic cross section. 

Fig. U shows the effect of hadronization on elliptic flow. Two extreme 
cases are presented: independent fragmentation and Lund fragmentation via 
400 long (Ay = 10) strings. In the case of independent fragmentation, the 
produced hadrons move in slightly different direction than the original parton 
(jet cone effect), resulting in a moderate reduction of elliptic flow. Lund frag- 
mentation produces all hadrons with almost complete azimuthal symmetry [], 
leading to an order of magnitude reduction at low p±. 



5 Conclusion 

The elliptic flow results from the parton cascade MPcSi strongly suggest that 
the observed saturation of elliptic flow at p± ~ 2 GeV can be quantitatively 

b The azimuthal symmetry is not perfect because the gluons kink the strings, however, the 
kinks are very small. 
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reproduced from parton kinetic theory, if the initial parton density is high 
enough. In addition, the sensitivity of elliptic flow on hadronization scheme 
will provide a strong constraint on possible-hadronization scenarios. A more 
detailed study of these issues is in progress!!! 
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